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Abstract 

Long range correlations multiplicity- multiplicity, ^-multiplicity and p 2 — "p\ are stud- 
ied in the percolating colour string picture under different assumptions of the dynamics 
of string interaction. It is found that the strength of these correlations is rather insen- 
sitive to these assumptions nor to the geometry of formed fused string clusters. Both 
multiplicity-multiplicity and p^-multiplicity correlations are found to scale and depend 
only on the string density, p^-multiplicity correlations, which are absent in the indepen- 
dent string picture, are found to be of the order of 10% for central heavy ion collisions 
and can serve as a clear signature of string fusion. In contrast p\ — "p\ correlations turned 
out to be inversely proportional to the number of strings and so very small for relaistic 
collisions. 

1 Introduction 

The standard description of particle production in soft high-energy strong interactions is in 
terms of colour strings stretched between the participant hadrons, whose decay generates the 
observed particle spectrum. The number of strings grows with energy and atomic number of 
participants and reaches thousands in heavy ion collision at RHIC and LHC energies. With 
their density high enough, one expects that the strings begin to overlap and interact. Some 
time ago a scenario of this interaction was introduced based on the percolation phase tran- 
sition and existence of colour strings of high colour [1]. An immediate consequence of the 
percolationg string scenario is damping of the multiplicity and rising of the average trans- 
verse momentum as compared to the independent string picture. These consequences are 
in agreement with the first data obtained at RHIC [2]. It is to be stresed that within the 
percolating string scenario one may imagine different versions of of the string interaction 
dynamics. Overlapping strings may not interact at all. Then the consequence of overlapping 
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will be just formation of different spots in the interaction area with higher colour content 
corresponding to the sum of colours of the overalpping strings. These spots will act as in- 
dependent particle emitters ("overlaps, or " ministrings" ) . The number of ministrings will 
generally be much greater than that of initial strings. So in this scenario percolation and 
fusion of strings actually leads to a proliferation of particle emitters. An opposite scenario 
with a lot of interaction of overlapping strings assumes that their colour becomes homoge- 
neously distributed over the formed cluster, the latter as a whole becoming an independent 
particle emitter. In this scenario the number of particle emitters is evidently smaller than 
the number of initial strings and may become unity if all strings form a single cluster. It is 
remarkable that these two scenarios actually lead to practically the same predictions for such 
global quantities as the multiplicity and average transverse momentum. So to distinguish 
between them one has to study more detailed information about particle spectra. Immediate 
candidates are the long range correlations between multiplicities or multiplicities and average 
transverse momentum. In this note we study these observables in both discussed scenarios 
of string interaction. Our results show that unfortunately these long range correlations are 
not very sensitive to the choice of string dynamics either. However there is one advantage 
of studying the long-range correlations involving the average transverse momenta: these cor- 
relations are absent in the independent string model. So there presence is a clear signature 
of string fusion and appearance of spots with higher colour density than on the average. 
As we shall see, the correlations between the transverse momentum and multiplicity are of 
especial value, since they depend only on the string density and not on the total number of 
strings. Long range correlations between transverse momenta, on the contrary, are inversely 
proportional to the total number of strings and so are practically absent in central heavy ion 
collisions. They can only be expected to be visible in highly peripheral collisions or collisions 
involving light nuclei. 

Some comments are to be made to explain the application of the colour string approach 
to correlations. First we shall study long-range correlations, which correspond to choosing 
the correlating observables from different rapidity windows separated by a reasonably large 
rapidity interval. This is to exclude short range correlations which are always present and 
appear both at the string decay stage (internal correlations of string particle production) 
and final hadronization stage (resonance decays). Second, in the study of average transverse 
momentum we assume the chosen rapidity windows to exclude the fragmentation regions, 
where the correlations between the average transverse momentum and multiplicities may 
largely be a result of summing the transverse momenta of the parent partons from the colliding 
hadrons. In fact the latter explanantion is the standard one for the observed correlations 
[3]. However it does not apply to the central region, in which particle production becomes 
insensitive to the parent parton momenta. So hopefully the correlations in the central region 
result only due to string interactions. Finally we have to remember that the unperturbative 
colour string picture in the present version is oriented only to soft particle spectra. So our 
region of average transverse momenta has to be limited by values of the order of 1-=- 2 GeV/c. 
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2 General formalism 

In this subsection we shall present our basic probabilistic expressions for the correlations 
valid for any scenario of string interaction. 

In a given event colour strings occupy certain regions in the transverse interaction plane 
and may overlap forming clusters. In any scenario this generates independent emitters, either 
overlaps or clusters as a whole. We enumerate these emitters by subindex a = 1,2,3, ...M, 
where M is the total number of emitters. 

An event consists of the emission of m a particles from emitter a. The total number of 
emitted particles is 

m e = J2 m a (1) 

a 

The pt distribution of particles emitted from emitter a will be given by a certain function 
w a (p). Note that it corresponds to the average 

Pi = J d 2 pp 2 w a (p) (2) 

whose concrete value depends of the scenario for the interaction. The effective pt distribution 
in a given event is then 

We{p) = m a w a{p) (3) 

m e a 

and the average p\ is 

Pe = —J2 m <*Pl ( 4 ) 

e a 

Now we pass to taking averages over many events. This averaging may be divided in two 
steps. First we average over various events occurring with the same string geometry, that is, 
with the fixed overlap and cluster structure (" configuration" ) . Then we have to average over 
all possible configurations. 

To average over events at a fixed configuration we have to know the probability p a (m) for 
a given emitter to produce m particles. This probabilty has to lead to the average number 
of particles 

m a = ^2mp a (m) (5) 

m 

which again is to be taken according to the chosen scenario of string interaction. We get 
averages in a given configuration: 

fhc= Y[pc t (. m a)m e = J2fh a (6) 

{m a } a a 

P 2 c= J2Y[pa(m a )^-Ym a p 2 a (7) 

{m a } « me « 

and also the probabilities in a given configuration to emit m particles 

(8) 

{m a } a 
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and a particle with the transverse momentum p 

{m a } oc ^ a m a 

To study correlations we have to know double distributions. We assume that emissions 
in the forward and backward rapidity windows are independent. As a result the double 
distribution in the numbers of particles mj? and vtlb emitted in the forward and backward 
rapidity windows in a given configuration factorizes 

p c (m F ,m B ) = p F (m F )p B (m B ) (10) 

where p^ B (m) are given by (8) with emitter probabilities p F,B (m) to produce particles in 
the forward (F) or backward (B) rapidity windows. Similarly the double distribution in 
the number of particvles m F in the forward rapidity window and the particle transverse 
momentum p B in the backward rapidity window in a given configuration is a product 

Pc{m F ,p B ) = p F {m F )w B (p B ) (11) 

Here w B {p) is given by (9) with p a —> p B . Finally the double distribution in momenta p F 
and pb of particles emitted in the forward and backward rapidity windows factorizes as 

w c {pf,Pb) = w B {p F )w B (p B ) (12) 

To average over different configurations we have to sum these expressions over all con- 
figurations with the weight P c , which is the probability of a given configuration. We denote 
this operation as < >. In particular we find various conditional probabilities of interest. The 
conditional probability to find tub particles in the backward rapidity window provided one 
sees mp particles in the forward rapidity window is 

, l w {pc(m F ,m B )) . . 

< P (m B )> mF = {pF[mF)) (13) 

The conditional probability to find a particle with transverse momentum ps in the backward 
rapidity window provided one sees m F particles in the forward rapidity window is 

/ ( w (pc(m F ,p B )) . . 

< W {jPB)> mF = {pF{mF)) (14) 

Finally the conditional probability to find a particle with momentum pb in the backward 
rapidity window provided one sees a particle with momentum p F in the forward rapidity 
window is 

. f . (w c {pf,Pb)) 

< "(PB > PF = {wF(pF)) (15) 

Taking the averages with these probabilities we find our basic formulas for correlations. 
The average multiplicity in the backward window at a given multiplicity in the forward 
window is 

. rn . _ ({m B )cPc{m F )) . . 

< mB >mF ~ (p B (m F )) (16) 
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where (ms) c and p F (mF) are given by (6) and (8) for the two rapidity windows. The average 
transverse momentum squared in some rapidity window at a given multiplicity in a different 
rapidity window is 

.2. ((PB 2 )cPc(m F )) 
PB mF= (p F (m F )) (17) 

where (pb 2 )c is given by (7) with p a — > p^. Finally the average transverse momentum 
squared in some rapidity window at a given momentum in a different rapidity window is 

__ {{PB 2 ) c w F ip F )) nR x 
(< (p F )> 

where iu^ is given by (9) with p a — > These formulas serve as a starting point for our 
study. 

Expressions for the correlations can be substantially simplified if the emission probabilities 
p a (m) have the Poisson form 

Pa(m) = P^ a (m) = e~ m ^ (19) 

Then taking average in a given configuration can be done analytically. 
In particular we find with (19) (see Appendix 1.) 

Pc(m) = Prn c (rn) = e-™^ (20) 

m! 

P 2 c = ^rT,™"Pl ( 21 ) 

five 

and . 

w c {p) = —^2m a w a (p) (22) 

m c a 

For particles emitted in the forward (backward) rapidity windows one has to substitute in 
these formulas fh —> ffiF,B 

With this simplification, calculation of correllations reduces to taking the averages only 
over configurations. 



3 Multiplicity correlations 

Forward-backward (FB) correlations between multiplicities has long been studied both ex- 
perimentally [4] and theoretically [5-7]. On the theoretical level they demonstrate that the 
number of elementary emitters fluctuates. On this ground they have served as a confirma- 
tion of the colour string picture, showing that the number of strings fluctuates, with its 
average growing with energy. For the FB multiplicity correlations to exist, the interaction 
between strings is irrelevant: the correlations are fully present when strings are independent. 
Moreover, as conjectured in [8], interaction and fusion of strings lead to damping of these 
correlations, since the effective number of strings diminishes. As we shall see this conjecture is 
confirmed by calculations using Eq. (16) in any of the two scenarios for the string interaction. 
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3.1 Ministrings as emitters 

As mentioned in the Introduction the two scenarios differ in the intensity of the interaction 
between the overlapping strings. The ministring scenario assumes that this interaction is 
in fact absent: the overlap regions retain their form and serve as independent emitters of 
particle with characteristics determined by the total colour accumulated in the overlap. In 
this scenario different emitters can be labeled by two numbers a = {ni} where i = 1,2,... 
enumerates overlaps of n strings. The basic quantity which characterizes the ministring 
scenario is the average multiplicity of overlap {ni} [4,5]: 

m ni = a Q y/n— (23) 

Here S n i and <7q are transverse areas of the overlap and initial string respectively. Using (25) 
we obtain the average multiplicity in a given configuration as 

m c = no Y] Vn— (24) 

where S n = J2i Sm is the total area in which n strings overlap and uq is the average multi- 
plicity for the initial string. The average multiplicities in the forward and backward rapidity 
windows are given by the same formula with uq ^ fi Q ' , the forward and backward multi- 
plicities for the initial string: 



- F,B F,B /—Sn / nc \ 

"V =/V Z^v^ - (25) 

n a 

Putting (25) into (20) and using (16) we reduce calculation of the FB multiplicity correla- 
tions to averaging over the overlap geometries of the known expression depending on S n . This 
averaging involves, first, different ways in which string overlap with their total N number 
fixed and, second, different values of N distributed around a certain average value < N > 
determined by the energy and atomic number of participants. At large < N > and small 
do as compared to the total interaction area S one expects a scaling behaviour: appropri- 
ately chosen characteristics of the correlations depend only on the dimensionless percolation 
parameter 

< N > a 

V = g (26) 

This property allows to study the correlations by Monte-carlo simulations with comparatively 
low values of < N > . Values of n and < N > for S-S and Pb-Pb central collisions at different 
energies obtained by Monte Carlo simulations based on the fusing colour string model [1] are 
shown in Table 1. 

In our numerical caluclations we studied the quantity 



{m B ) 



= - 1 (27) 



as a function of 

* = ^ (28) 
(m F ) 
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In our simulations we used both the homogeneous distribution of strings in the transverse 
interaction area and the inhomogeneous one, which follows the nuclear profile function T(b). 
No significant qualitative difference was found between these two choices. For simplicity we 
present here our results for the homogeneous distribution. We chose fi ' =1. We used the 
Poisson distribution in the number of strings N with < N >= 25 and 50. We have checked 
that the results are practically independent of < N > in the region 0.3 < x < 3, so that 
F^n(x) exhibits the expected scaling behaviour for physically relevant values of mj?. 

Our results for values of r] = 0.5 -j- 3 show that in the whole interval < x < 3 
is a monotonously rising function of x, rather close to linear and crossing zero at x = 1. Its 
slope at x = 1 (and < N >= 50) (the "correlation coefficient") 



dx 



is shown in Fig. 1 by a solid line. As we observe it steadily falls with 77, in full accordnace 
with the idea that the correlations diminish as a consequence of string fusion. At very high 
rj the slope seems to flatten against the value around 0.1. At 77 — ► the slope tends to the 
value 0.5, which corresponds to the independent string model with the assumed value of fj,o 
(see Appendix 1). Note that 77 ~ 1.12 corresponds to the critical value for the percolation 
phase transition. However our curves do not show any peculiarity at such values of 77. 

In Fig. 2 we compare b^—^ for the homogeneous distribution of strings and the inhomo- 



geneous one which follows the nuclear ptofile function T(b) oc y R\ — b 2 (corresponding to a 
constant nuclear density). As mentioned, the difference is of no practical importance. 

3.2 Clusters as emitters 

An alternative scenario with a strong interaction between overlapping strings assumes that the 
colour becomes homogeneously distributed over the whole cluster formed by the overlapping 
strings. The emitters are now clusters themselves and index a now labels just different 
clusters. The average multiplicity of a cluster is chosen to adjust to the colour distribution 
[9]: 

fh a = fio\- £L - £L (29) 
V °o 

where n a and S a are the number of strings and area of cluster a. For the forward and 
backward rapidity window (29) transforms into 

(30) 

Eqs. (29) and (30) express the physical contents of the cluster scenario for multiplicity corre- 
lations and their difference from Eqs. (24) and (25) shows the effect of the strong interaction 
between overlapping strings. The remaining procedure to calculate the correlations does not 
change as compared to the ministring case. 

We calculated the same function F^_ M (x) (Eq. (29)) with (32) by Monte Carlo simula- 
tions. Comparison of our results for < N >= 25 and 50 again showed a satisfactory scaling 
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in < N > for x = 0.3 3 and a weak dependence on the way the strings are distributed 
in the interaction area. In spite of a very different dynamics, we have found no appreciable 
difference with the former case ov overlaps as emitters. The slopes of at x = 1 

are shown in Fig. 1 by a long-dashed line (again for < N >= 50). They are systematically 
higher than for overlaps but the difference does not exceed 10%, which is of no practical 
importance, taking into account simplifications involved at the basis of the model. 

3.3 Fixed number of emitters 

Technically taking the average over all geometrical distributions of strings in the transverse 
area is a formidable task, which can be realistically achieved by Monte Carlo simulations 
for a reasonable time only for a limited number of strings, substantially smaller than this 
number in heavy-ion collisions at RHIC and especially LHC energies. On the other hand, it 
is well-known that different ways of averaging often lead to the same or practically the same 
averages, since the physical picture of fluctuating variables is often basically the same. This 
motivates searching for a simplified picture of string fusion, which, on the one hand, is not 
very different from the discussed above by its physical implications and, on the other hand, 
avoids geometrical averaging. In this subsection we propose such a picture. 

Physically the basis of string fusion model consists of appearance of various emitters 
homogeneously distributed in the transverse plane with the number of overlapping strings 
varying for 1 to JV with a certain probability. One can model this physical situation by 
assuming a fixed number M of emitters "cells"), each one of them (emitter a) equivalent 
to a certain number n a of completely overlapped strings including n a = 0, in which case 
there is no emission at all. So effectively one also has a varying number of emitters with 
different colours. Averaging over configurations is achieved by distributing each n a around 
some average value n a with a certain probability P a (n a ). The homogeneous distribution of 
strings in the transverse plane corresponds to equal averages: n a = n . One can model an 
inhomogeneous distribution of strings by distributing the M emitters in the transverse plane 
and choosing n a in accordance with the nuclear profile density. 

In this model ("cell model" [10]) we have for the average number of particles from cell a 



If we assume that the cell probabilities p a (m) are Poissonian then correlations can be calcu- 
lated using our formulas of Section 2 in exactly the same way as before. For the final averaging 
over configurations (that is over all n a , a = 1, ...M) we have taken the probabilities P a (n a ) 
to be also Poissonian with a common average value n = rj to correspond to the previous 
geometrical picture. Note that in this model the average value of strings is < N >= Mr). 

The results of calculations of F^^ix) in this model for < N >= 50 and various r] show 
that they are very close to both previous models, especially at relatively large rj when they 
are nearly identical to those in the ministring picture. The slopes of F^-^(x) at x = 1 for the 
cell model are shown in Fig. 1 with a short-dashed line. So as expected, this simplified model 
imitates the geometrical model of ministrings almost ideally. With that, it is much simpler: 
its realization by Monte Carlo simulations requires computer time more than an order of 
magnitude smaller than for ministrings with equal < N >. Apart from this, the model 




(31) 
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admits some analytic estimates for small or large values of 77, which reveal basic properties 
of the correlaltions in these asymptotical regions (see Appendix 2.). 



4 Correlations between average transverse momentum and 
multiplicities 

According to Eq. (17) to calculate the correlations < p\ > — fi we have to know the averages 
p 2 c in the backward window in a given configuration. For the Poissonian emitters they can 
be calculated via Eq. (21) if one knows these averages for a given emitter a (Eq. (2)). The 
concrete values of p 2 a depend on the scenario of string interaction. They were analyzed in 
[8,9] wherefrom we borrow their expressions presented below. 

Both for the ministring and cell scenarios the value of p 2 a for the individual emitter are 
determined only by the the number of strings in the overlap or cell. For ministrings 

V 2 m = P 2 oV^ (32) 

and for the cell model 

P 2 a =plVn^ (33) 

For the cluster scenario p 2 a depends on the color density inside the cluster and thus also 
depends on its area S a : 

2 2 l n a&0 , QA y. 

Pa=Po \~S^ ^ ' 

In these formulas p\ is the average transverse momentum squared for the initial string. 

With these expressions for p 2 a we calculated (p 2 B )m F according to Eq. (17) by Monte Carlo 
simulations for the three scenarios considered: ministrings, clusters and cells. We studied 

F p >-» - - 1 (35) 

considered as a function of the same variable x, Eq.(28). As before we took < N >= 25 and 
50 and 77 in the region 0.5 -j- 3. We found that F p 2_ )l (x) again exhibits a satisfactory scaling 
property: it results practically independent of < N >. As for /j — /x correlations F p 2_ fl (x) is 
found to be a nearly linearly rising function of x, crossing zero at x = 1. However, contrary 
to n — /i correlations, the slopes for p 2 — /j correlations rise with 77. Slopes of F p 2_ fl (x) at 
x = 1 

/ dFi_ Jx)\ 



dx j 

/ x=l 



are shown in Fig. 3 for the three studied scenarios with the homogeneous distribution of 
strings. One observes that the overlap and cluster scenarios again lead to practically identi- 
cally results. The cell model imitates these two physical scenarios quite satisfactory at large 
r] > 2 4- 2.5. At smaller r\ its slopes lie 30-^40% below. At large r\ slopes in all three models 
seem to converge around 0.1. It is remarkable that this value coinsides with the corresponding 
limit for the slopes in /i — /j correlatioms (cf. Fig. 1). At 77 — ► all slopes vanish, which 
corresponds to absence of p 2 — fi correlations for independent strings. 
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In Fig. 4 we compare the slopes in the overlap scenario for the homogeneous distribution 
of strings and the same realsitic one, as in Fig. 2. Again the difference results of little 
importance. 

It has to be stressed that absolute values of the slopes are not too small, of the order 
of 0.1 in a wide range of rj. Such slopes can hopefully be measured experimentally. There 
existence would be a clear signature of string fusion. 



*«(p) = v :z ' ( ) (36) 



5 Pb — Pf correlations 

The last type of correlations we shall study are those between transverse momenta of particles 
emitted in the forward and backward rapidity windows. They can be measured by the average 
transverse momenta squared of particles in the backward rapidity window at an observed 
transverse momentum in the forward rapidity window, Eq. (18). To calculate this in our 
colour string model we have to know the particle spectrum w a (p) for an individual emitter. 
This spectrum has to be chosen to lead to the averages (2), which have been discussed for 
each of the models for string interaction in the previous section. We choose the distribution 
w a (p) in the form prompted by the standard fit to the experimental spectrum for the proton 
target [11] 

(K -l)(«-2) / Pa 

2ttpI \p + p a/ 

Here p a is given by Eqs. (32)-(34) for the three models of string interaction considered, 
Pa = y/pj, 

k =19.7 -0.86 In E, = (GeV /c) 2 (37) 

and E is the cm. energy in GeV. 

With this distribution we calculated (18) for the three considered scenarios averaging 
over configurations by Monte Carlo simulations with different average number of strings 
< N >= 25 and 50. Calculations show that for these correlations the scaling function 
independent of < N > turns out to be 

F p _ p =<N> - P \ >PF - 1 (38) 
< p\ > 



As argument it is convenient to choose dimensionless 

PF 



\I<P 2 F> 



(39) 



We have calculated F p - p (y) for the for different r\. In all cases F p - p (y) is found to be a 
linearly rising function of y crossing zero at y = 1. The slopes of F p - P (y) at y = 1 



bp—p 



dFp~ P (y) 

dy 
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are shown in Fig. 5 as functions of rj for the three studied models. They rise with rj similarly 
to the slopes in p 2 — jj, correlations. Again overlaps and clusters lead to practically identically 
results. Their imitation by cells gives rather satisfactory results, especially at rj > 2 where 
the cell slopes lie higher than those from the other two scenarios by some 10%. At lower rj 
this difference rises to 30-^40 %. 

Note that the presence of the factor < N > in (38) together with the calculated values of 
the scaling function show that for realistic collisions with high enough values of rj the magni- 
tude of the correlations is extremely small, which excludes their experimental observation. As 
follows from our results the correlations can be observed in practice up to < N >= 10 15, 
which is unrealistic for rj > 1. 

6 Conclusions 

This work was primarily motivated by search of appropriate observables which could shed light 
on the details of dynamics of colour string fusion conjectured to occur in heavy-ion collisions. 
The results of our study show that pair correlations are unfortunalely rather insensitive to 
these details. Moreover they demonstrate that even the geometric picture of string fusion and 
percolation is not very essential for the found correlations. A very simplified picture, which 
only takes into account existence of many emitters with different colours, leads to results not 
very different from the realistic geometrical overlapping of colour strings. So the correlations 
basically reflect only the possibility of spots with higher colour field in the colliding nuclei 
and strongly depend only on their spatial density characterized by the parameter rj. 

From a certain point of view this is a positive aspect, since predictions for the correlations 
do not vary with specific assumptions made about the string fusion dynamics. The behaviour 
of the correlations with r\ emerges as a clear signature for the basic new phenomenon char- 
acteristic for string fusion: formation of emitters with higher colour density. Both ji — ji 
and p 2 — jx correlations look quite promising from this point of view. In particular p 2 — ji 
correlations arise totally due to string fusion in our picture; they are absent with independent 
strings. However, as mentioned in the Introduction, this is only true for particles produced 
strictly in the central region, since in the fragmentation regions obvious p 2 — ji correlations 
are caused by the accumulation of the string ends momenta. So to see the effect of string 
fusion one should choose both rapidity windows as far as possible from the rapidity limits. 
As to p — p correlations, our calculations have shown that they diminish with the number 
of strings and appear to be negligible for realistic heavy-ion collisions except at the highly 
peripheral collisions. 

Our calculations have been made for a homogeneous distribution of strings in the trans- 
verse interaction plane. Calculations for a realistic distribution, taking into account the 
varying string density as a function of impact parameter, are now in progress. However the 
found insensitiveness of our results to the details of geometry allows to expect that the change 
introduced by a realitic geometry will be not very essential. 
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8 Appendix 1. Averages with the Poisson distribution 

We start with p c (m) determined by (8). Presenting the Kronecker symbol as a contour 
integral around the origin we find 

Pc(m) = £ / ^- z z ~ m+Eama I[pa( m «) ( 4 °) 

{m a } " 

We have 

£ z m «P{m a ) = e^" 1 )*- (41) 

rn a 

so that 

M™) = /^<*- 1, *=e- ffl 'S W2) 

So we find (20). 

Now take average (7). We present it as 

P 2 c = £ — £ II P«( m «) 5 m,^ Q m Q £ ( 43 ) 

The internal sum over all m a at their sum fixed can be represented in the same way as (40): 

x = £ / ^r~ m+E - ma n^KOE^ (44) 

{m a } a a 

Take some particular term in the last sum, say a = 1 Then sums over all m a , a = 2, 3, ...M 
will give the same factor (41). The sum over mi will however contain an extra factor mi: 

E z mi miP(mi) = e^™ 1 = zrhie^ 1 ^ (45) 

mi 

So we find 

X = E m aP le-^ f ^e* = ^p c (m) E (46) 

where p c (rn) is given by (20). Putting this into (43) we find that the factor m in (46) cancels 
the denominator in (43). Using then 

EpM = 1 

m 

we find (21). 

Derivation of (22) is done in exactly the same way. 
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9 Appendix 2. Analytic estimates in the cell scenario 

As mentioned, the cell scenario admits explicit analytic estimates [10] valid in the limit of 
large or small values of r). 

Estimates at large r\ are based on the asymptotic equivalence of the discrete Poisson 
distribution and continuous Gaussian distribution 

P n {n) ~ G n (n) = -L= e ~^ (47) 

valid in the limit n » 1. 

Let us first study the simplest case of independent strings (without fusion). Then Eqs. 
(31) and (33) are changed to 

m a = ^n a , vl=vl (48) 

and similarly for averages in the two rapidity windows. To simplify in this Appendix we shall 
take \i F = /j,q and denote this common multiplicity by fio- From the second of Eqs. (48) 
it follows that there will be no p 2 — ]x correlations as expected. The configuration averages 
fh F,B = fio J2 a n a = Mo-^c wm depend only on the total number of strings N c in a given 
configuration. Averaging over configuration will be reduced to averaging over the number 
of strings in different configurations. If the number of strings in each cell is distributed 
according to the Poisson distribution with average i] then the total number of strings will 
also be distributed according to the Poisson distribition with the average Mr] (see Appendix 
1). So at fixed r] the overall average number of strings is related to M as 

(N) = Mr] (49) 

Thus, substituting the two Poisson distributions by Gaussians we find from (16) 



' dN N-We-W'"*) 



{m B )m F = ° r oo JAr „-1/?.-MW m ^ ( 50 ) 



where 

^ mF) = 2(N) + 2^ N (51) 
We estimate the two integrals by the saddle point method. Then we find 

{m B )m F = fJ, No(m F ) (52) 
where iVo is the solution of the equation d<j>/dN = 0. In terms of scaling quantities 

X= (^f)' z = Jn) (53) 

the equation for the saddle point takes the form 

z 3 -z 2 = ^]i (x 2 - z 2 ) (54) 
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Solving this equation for zq{x) one can find function (mB) mp (mj?) at all mp. The correlation 
coefficient 6^-^ can however be determined explicitly from (54): 

(55) 

\dxj x=1 Mo + 1 

Taking /io = 1 we have 6^=^ = 0.5 in full agreement with the limit of curves shown in Fig. 1 
at 7] — > 0, which corresponds to independent strings. 

Now we pass to a physically more interesting case when strings fuse, which in the cell 
scenario corresponds to Eqs. (31) and (33). Introducing a quantity for a given configuration 

r c = V™a ( 56 ) 

a 

we find 

fh c = fi r, p 2 c = P l-^- (57) 

so that 

{mB)mF ~ M ° (P,Am f )) (58) 



and 



The average is taken over the product of Poisson distributions in n a each with with the 
average rj. 

We consider r\ » 1 and substitute all Poisson distributions with Gaussian ones to get 

f°°T7 tin r r ^ X l 2 p -4>{n a ,m F ) 
I \ JO lla u ' i a'c'c c . . 

(m B ) mF = Mo^— — _ 1/2 ~ (60) 



and 

2 2 Jo 00 Ua dn a (N c /r c )r- 1/2 e-^ m ^ 

\PB/m F — Po — T— " ~ : 

io EL dn a r c ' e-ttn*,™*) 

where 

, / n (n a - r]) 2 (m F - por c ) 2 

^^) = E^r- + 2/iore 

In the saddle point approximation we find 



(61) 
(62) 



(m B )m F = Por , (p 2 B)m F =pl— (63) 

^0 

where N an d r are values of N c and r c at the saddle point defined by the equations 

d<P(n m F ) = a = 
dn a 



15 



Introducing scaled variables 
and a constant involving r\ 



{m F Y Za ~ V r, ' (b5j 



Mo 



4y^ 

we can write the saddle point equations in the form 



(66) 



zl-z a = a(x 2 ^P--l\ , a = l,...M (67) 



In terms of z a here r c = ^/rjJ2a z a- Due to symmetry in a obviously for the solution z a = z 
where z satisfies a single equation 

, 3 -z = a(jJ-l) (68) 

This equation defines z = z(x), in terms of which one finds 

F^(x) = F p 2^ = z{x) - 1 (69) 

So at large values of r\ both ft — fi and p 2 — \x correlations are described by the same scaling 
function of x. 

At x = 1 Eq. (68) possesses an obvious trivial solution: 

x = 1, 2=1, r = Mjfj, //N = (N) = fi M (70) 

This is sufficient to calculate the (identical) correlation coeficients for both fi — fi and p 2 — ft 
correlations: 

W = V - = ^TI^ (71) 

It slowly falls with r] Of course one should have in mind that thexe are only asymptotic 
estimates, valid at sufficiemtly high r]. To see their validity region we present in Fig. 6 the 
slopes b in fi — fi and p 2 — fi correlations as functions of 77 with fiQ = 1 calculated by Monte- 
Carlo simulations, together with the common asymptotic curve (71). As one observes the 
asymptotic curve gets more or less close to the exact ones starting from rj > 4 -j- 5. These are 
very high values from the physical point of view, attainable only at energies in the multi-TeV 
region for the heaviest nuclei and central collisions. 
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11 Table 

Values of rj and average number of strings < N > for central S-S and Pb-Pb 
collisions at different energies 

S-S scattering (b = 0) 





V 


< TV > 


19.4 


0.40 


126 


62.5 


0.52 


164 


200 


0.65 


204 


546 


0.80 


252 


1800 


0.97 


308 


7000 


1.23 


390 



Pb-Pb scattering (b = 0) 





V 


N 


19.4 


1.08 


1190 


62.5 


1.39 


1536 


200.0 


1.60 


1800 


546.0 


2.03 


2240 


1800.0 


2.46 


2720 


7000.0 


3.03 


3350 
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Figure 1: Correlation coefficient b for fi — fi correlations in the three scenarios studied. 




Figure 2: Correlation coefficient b for n — \i correlations in the overlap scenario for the 
homogeneous distribution of strings and the one proportional to T(b) oc J R\ — b 2 . 




Figure 4: Correlation coefficient b for p 2 — \i correlations in the overlap scenario for the 
homogeneous distribution of strings and the one proportional to T(b) oc J R\ — b 2 . 




Figure 6: Correlation coefficients for /x — \i and p 2 — jjl correlations in the cell model together 
with the asymptotics (71). 



